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Radial Profiles of Gas in Late-Type Disk Galaxies 

Curtis Struck-Marcell, Iowa State University 

x, I. Review of Published Data 


The azimuthally averaged HI distribution, and the total gas density distribution derived from HI and 
CO observations (and assuming I% 2 = 2.8 x 10 20 I c0 ) as a function of radius in several nearby, 
early-type disks are examined in this study. 

A. Z°e 1/r where the Rotation Curve is Flat 


First pointed out by Quirk (1972), and recently discussed by Kennicutt (1989) this conclusion is 
evident in the gas distributions considered here, and in the data of Wevers' (1984) thesis. For 
example, the figure shows the deviation of the azimuthally averaged HI surface density profiles from a 
best-fit 1/r profile in Wevers' galaxies. Only data from positions with circular velocities within 5% of 
the average value in the flat part of the rotation curve are plotted, only those galaxies with more than 3 
such points were included. The curve for NGC 628 also includes the molecular gas observed by Adler 
and Lizst (1989). For the other sample galaxies the molecular gas either lies interior to the flat part of 
the rotation curve or a CO distribution has not been published. It is apparent that most of the data fall 
within about 25% of the mean 1/r profile, except in the outermost parts of some galaxies. The larger 
surface density deviations in these outer regions frequently correlate with deviations from flatness in the 
rotation curve. 


B. The Gas Distribution in the Inner Disk is More Complex 

The available data for the inner regions of the disks are more sparse and their interpretation is more 
problematic. We can only make a few tentative generalizations: 

1. When the central gas surface density is high (e.g. X > 30 solar masses/pX) the 1/r profile continues 
from the outer regions to deep into the central regions (e.g. IC 342, NGC 5055, NGC 6946). 

2. In the other galaxies there is still a central peak in X(r). In the cases considered here, this peak is 
clearest in M33 and NGC 2403, where it lies within the region with v^ < 60 km/s. 

3. If the central peak is well interior to the X 1/r region, then there is a flat, X «= constant, 
transition region (M33, NGC 2403). In NGC 628 the center may be primarily of this flat form, since 
the central rise is very modest. 

n. Theoretical Considerations and Constraints 


A. Azimuthally Averaged Hydrostatic Equilibrium (for a continuum gas) 
We have accretion disk equations like those in the review of Pringle (1981). 


1. Mass Continuity, implies 9X/9t = 0, radial velocity v r = 0, and self-diffusion term 9 2 (rX)/9r 2 = 0, 
which yields X = a/r + b. with a,b = constant. 

2. Equation of State: Isothermal cloud gas, with sound speed c = constant, may be an adequate 
approximation. 

3. Momentum Transport: 

a) Radial Force Balance (in a potential dominated by collisionless matter). The azimuthally averaged 
pressure gradient force is 


= Xo-ic 2 ) 


rX 


dr 


rX 


dr 


Case i) X« 1/r. The pressure gradient force is zero, and centrifugal force balances gravity, just as 
for the stars. 
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Case ii) X = constant. The averaged pressure gradient force goes as 1/r. In the flat rotation curve 
region, v e = constant » c, and both gravity and centrifugal force also go as 1/r, so pressure is 
negligible throughout the region. In a solid body region, v 0 r, so as r decreases the pressure 
increases relative to the centrifugal force. 


b) Azimuthal Force Balance 


The only force is viscosity, so in equilibrium it must be zero, then, 


vXr 3 -^ 

dr 


— = constant, 


which gives a relation between X and v 0 . E.g., 

if v e = constant, then vXr = constant (or X 1/r for v = constant). 

if v 0 r, then d(vo/r)/dr = 0, and X is unrestricted. 

if X = constant, then either v 0 « r or v 0 1/r (unphysical). 


The role of (spiral) waves in transporting angular momentum, and thus influencing radial structure is 
not clear, and will not be considered here. 


B. Local Gravitational Stability (Toomre Q) 

Larson (1988) and Kennicutt (1989) have recently summarized the theoretical indications and 
observational evidence, respectively, showing that Toomre's (1964) gravitational stability parameter Q 
= ck/tiGX is of order unity throughout the star-forming part of a galaxy disk. In the isothermal, v 0 = 
constant part of a disk this also implies X « 1/r, consistent with the other constraints, but also setting 
the value of X. 


C. Secular Kinetic Equilibrium 

Large clouds may have long mean free paths (k ~ 1 kpc.) in galaxy disks, so kinetic (random walk) 
effects must be considered. Diffusive mass fluxes between adjacent annuli of width X. must balance, 
soXrc = constant. (For practical purposes this is equivalent to zero self-diffusion.) Momentum 
Flux balance between adjacent annuli implies, Xrc v 0 = constant. This condition is not equivalent to 
minimizing hydrodynamic viscous shear. E. g. the latter is zero when v 0 «= r, and X «= 1/r, but the 
present condition is not satisfied! 


III. Conclusions 

- Both the available data and theoretical considerations suggest that the unique equilibrium surface 
density profile in the flat-rotation-curve region of the disk is X « l/r. If this profile is also stable to 
non-axisymmetric disturbances, then this region may be regarded as a gas reservoir for the galaxy. 

- The existence of equilibrium or quasi-steady profiles in the inner disk, where v 0 increases with r is 
more problematic. If v e r in this region, then there are two hydrodynamic equilibrium forms: X « 
1/r or X = constant. There are examples of both forms in the observations. However, neither is a 
kinetic equilibrium form, so we can expect diffusive processes to funnel material from the reservoir 
region inwards. At the same time, the critical density for local gravitational instability is constant in 
regions where v 0 r, so if X « X cr jj, star formation processes may oppose the kinetic ones. We do 
not yet understand the outcome of this interplay of forces. 
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Figure la-c. Shows the deviation of the azimuthally averaged gas surface density profiles from 
best-fit 1/r profile using the data from Wevers (1984) as described in the text. 
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